
********************************************************************************
********************************************************************************
*** Replication materils for: "Economic Sanctions and Food Consumption: Evidence from Iranian Households" ///
*** forthcoming in International Studies Quarterly
**** Babak RezaeeDaryakenari, Peyman Asadzade, and Cameron G. Thies
*** Part 2: regression analysis ************************************************
*********************************************************************************


** set the workign directory 

cd "~"

import delimited "Data_cl_Regressions.csv", clear

gen Sanction_OilComprehensive=0

label variable Sanction_OilComprehensive "Comprehensive sanction" 
** 2010 US EU sanction until the Irandeal become effective
*** We lagged the sanction variable as it takes sometime to come to effect!
replace Sanction_OilComprehensive=1 if year>=1392 & year<=1395
replace Sanction_OilComprehensive=1 if year>=1398

*Log income variable
gen income_total1000k_log=log10(1+income_total1000k)

* keep the years after 1370 (persian calendar), when is the focus of the analysis, which the years after the end of Iran-Iraq war 
keep if year>=1370

gen pre_period=.

replace pre_period=0 if year<=1372
replace pre_period=1 if year<=1376 & year>1372
replace pre_period=2 if year<=1380 & year>1376
replace pre_period=3 if year<=1384 & year>1380
replace pre_period=4 if year<=1388 & year>1384
replace pre_period=5 if year<=1392 & year>1388
replace pre_period=6 if year<=1396 & year>1392
replace pre_period=7 if year<=1400 & year>1396


***
* Lable variables

label variable qty_kg1_hh   "Flour" 
label variable qty_kg2_hh   "Grain" 
label variable qty_kg3_hh   "Bread" 
label variable fruitvegtable   "Fruits \& vegetables" 
label variable qty_kg5_hh   "Red meat" 
label variable qty_kg6_hh   "Poultry" 
label variable qty_kg8_hh   "Seafood"
label variable qty_kg8_hh   "Dairies \& eggs" 

**********

*** Summary table
set matsize 11000
outreg2 using "Tables\Table_A1_SummaryTable.rtf"  ///
            , replace sum(log) label ///
			keep(qty_kg1_hh qty_kg2_hh qty_kg3_hh fruitvegtable qty_kg5_hh qty_kg6_hh qty_kg8_hh milkeggdairies ///
			Sanction_OilComprehensive hhemployed hhmale hhfemale hhagemean income_total1000k_log urban quartile)

outreg2 using "Tables\Table_A1_SummaryTable.tex"  ///
            , replace sum(log) label ///
			keep(qty_kg1_hh qty_kg2_hh qty_kg3_hh fruitvegtable qty_kg5_hh qty_kg6_hh qty_kg8_hh milkeggdairies ///
			Sanction_OilComprehensive hhemployed hhmale hhfemale hhagemean income_total1000k_log urban quartile)

***************************
			
eststo clear
** reghdfe method is faster and also cluster the standard errors, see http://scorreia.com/software/reghdfe/

eststo: reghdfe qty_kg1_hh i.Sanction_OilComprehensive  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg1_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Flour", size(3)) saving(mrg_qty_kg1_hh.gph, replace)


eststo: reghdfe qty_kg2_hh i.Sanction_OilComprehensive hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg2_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Grain", size(3)) saving(mrg_qty_kg2_hh.gph, replace)


eststo: reghdfe qty_kg3_hh i.Sanction_OilComprehensive hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg3_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Bread", size(3)) saving(mrg_qty_kg3_hh.gph, replace)


eststo: reghdfe fruitvegtable i.Sanction_OilComprehensive hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store fruitvegtable
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Fruits & vegtables", size(3)) saving(mrg_fruitvegtable.gph, replace)


eststo: reghdfe qty_kg5_hh i.Sanction_OilComprehensive  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg5_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Red meat", size(3)) saving(mrg_qty_kg5_hh.gph, replace)


eststo: reghdfe qty_kg6_hh i.Sanction_OilComprehensive  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg6_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Poultry", size(3)) saving(mrg_qty_kg6_hh.gph, replace)


eststo: reghdfe qty_kg8_hh i.Sanction_OilComprehensive  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg8_hh
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Sea food", size(3)) saving(mrg_qty_kg8_hh.gph, replace)


eststo: reghdfe milkeggdairies i.Sanction_OilComprehensive  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store milkeggdairies
margins i.Sanction_OilComprehensive
mplotoffset, scheme(plotplain)  offset(0.02) xtitle("") ytitle("") title("Dairies & eggs", size(3)) saving(mrg_milkeggdairies.gph, replace)


esttab using "Tables\Table1_Tbl_Main.rtf", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE , fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted  ///
		/// drop(*.province) ///
		 nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001) ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)	
					
esttab using "Tables\Table1_Tbl_Main.tex", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE, fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted  ///
		/// drop(*.province) ///
		nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001)  ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)						
					

gr combine mrg_qty_kg1_hh.gph mrg_qty_kg2_hh.gph mrg_qty_kg3_hh.gph mrg_fruitvegtable.gph mrg_qty_kg5_hh.gph ///
		   mrg_qty_kg6_hh.gph mrg_qty_kg8_hh.gph mrg_milkeggdairies.gph, ///
		   commonscheme scheme(white_piyg) plotregion(fcolor(white)) graphregion(fcolor(white)) ///
		   l1(Household Consumption per Person per Month) b1(Comprehensive Sanction) /// 
		   col(4) iscale(.8) xsize(12) ysize(7) imargin(-4 -4 -4 -4) 
		   
gr export "Figures\Figure3_AllMargins.eps", replace 
gr export "Figures\Figure3_AllMargins.pdf", replace 
		   
************ Regression results in subgroups

** By urban rural

eststo clear

eststo: reghdfe qty_kg1_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg1_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Flour", size(3)) legend(order(3 "Rural" 4 "Urban") rows(1)) ///
  saving(mrg_qty_kg1_hh.gph, replace)


eststo: reghdfe qty_kg2_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg2_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Grain", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_qty_kg2_hh.gph, replace)


eststo: reghdfe qty_kg3_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg3_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Bread", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_qty_kg3_hh.gph, replace)



eststo: reghdfe fruitvegtable i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store fruitvegtable_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Fruits & vegtables", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_fruitvegtable.gph, replace)



eststo: reghdfe qty_kg5_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log  ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg5_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Red meat", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_qty_kg5_hh.gph, replace)



eststo: reghdfe qty_kg6_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg6_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Poultry", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_qty_kg6_hh.gph, replace)


eststo: reghdfe qty_kg8_hh i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg8_hh_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Sea food", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_qty_kg8_hh.gph, replace)


eststo: reghdfe milkeggdairies i.Sanction_OilComprehensive##i.urban hhemployed hhmale hhfemale hhagemean income_total1000k_log ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store milkeggdairies_RU
margins Sanction_OilComprehensive#urban
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) ///
  offset(0.02) xtitle("") ytitle("") title("Dairies & eggs", size(3)) legend(order(3 "Rural" 4 "Urban")) ///
  saving(mrg_milkeggdairies.gph, replace)


esttab using "Tables\TableA2_Tbl_UrbanRural.rtf", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE, fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted  ///
		/// drop(*.province) ///
		 nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001)  ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)	
					
esttab using "Tables\TableA2_Tbl_UrbanRural.tex", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE, fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted ///
		/// drop(*.province) ///
		 nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001)  ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)						

grc1leg2   mrg_qty_kg1_hh.gph mrg_qty_kg2_hh.gph mrg_qty_kg3_hh.gph mrg_fruitvegtable.gph mrg_qty_kg5_hh.gph ///
		   mrg_qty_kg6_hh.gph mrg_qty_kg8_hh.gph mrg_milkeggdairies.gph, ///
		   commonscheme scheme(white_piyg) plotregion(fcolor(white)) graphregion(fcolor(white)) ///
		   l1(Household Consumption per Person per Month) b1(Comprehensive Sanction) /// 
		   col(4) iscale(.8) xsize(12) ysize(7) imargin(-4 -4 -4 -4) 

gr export "Figures\Figure4_AllMargins_RU.eps", replace 
gr export "Figures\Figure4_AllMargins_RU.pdf", replace 
		
***** By income quartile

eststo clear

eststo: reghdfe qty_kg1_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg1_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Flour", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg1_hh.gph, replace)


eststo: reghdfe qty_kg2_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg2_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Grain", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg2_hh.gph, replace)


	
eststo: reghdfe qty_kg3_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg3_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Bread", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg3_hh.gph, replace)



eststo: reghdfe fruitvegtable i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store fruitvegtable_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Fruit & vegtable", size(3)) legend( rows(1)) ///
  saving(mrg_fruitvegtable.gph, replace)


	
eststo: reghdfe qty_kg5_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg5_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Red meat", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg5_hh.gph, replace)


	
eststo: reghdfe qty_kg6_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg6_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Poultry", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg6_hh.gph, replace)



eststo: reghdfe qty_kg8_hh i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store qty_kg8_hh_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Sea food", size(3)) legend( rows(1)) ///
  saving(mrg_qty_kg8_hh.gph, replace)



eststo: reghdfe milkeggdairies i.Sanction_OilComprehensive##i.quartile  hhemployed hhmale hhfemale hhagemean income_total1000k_log i.urban ,  absorb(admin_2code pre_period) vce(r)
estadd local locationFE "Yes"
estadd local timeFE "Yes"
estimates store milkeggdairies_Q
margins Sanction_OilComprehensive#quartile
mplotoffset, scheme(white_piyg) plot1opts(lcolor(gs14) msymbol(D)) plot2opts(lcolor(gs14) msymbol(Oh)) plot3opts(lcolor(gs14) msymbol(S)) plot4opts(lcolor(gs14) msymbol(Th)) ///
  offset(0.02) xtitle("") ytitle("") title("Dairies & egg", size(3)) legend( rows(1)) ///
  saving(mrg_milkeggdairies.gph, replace)



esttab using "Tables\TableA3_Tbl_Q.rtf", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE, fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted ///
		/// drop(*.province) ///
		 nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001)  ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)						

esttab using "Tables\TableA3_Tbl_Q.tex", replace se varwidth(25) level(95) label ///
		title(Table 1: Estimation results) ///
		s(N ll bic locationFE timeFE, fmt(%10.2f) ///
		label("N" "Log-likelihood" "BIC" "Location FE" "Presidential period FE")) ///
		nobaselevels interaction("X") noomitted ///
		/// drop(*.province) ///
		 nonumbers  mtitles("Flour" "Grain" "Bread" "Fruits \& vegtables" "Red meat" "Poultry" "Sea food" "Dairies \& eggs") ///
		compress nogap	star(* 0.05 ** 0.01 *** 0.001)  ///
		b(%5.2f)  ///
		coeflabels(_cons "Intercept"  ///
					)		

grc1leg2   mrg_qty_kg1_hh.gph mrg_qty_kg2_hh.gph mrg_qty_kg3_hh.gph mrg_fruitvegtable.gph mrg_qty_kg5_hh.gph ///
		   mrg_qty_kg6_hh.gph mrg_qty_kg8_hh.gph mrg_milkeggdairies.gph, ///
		   commonscheme scheme(white_piyg) plotregion(fcolor(white)) graphregion(fcolor(white)) ///
		   l1(Household Consumption per Person per Month, size(3)) b1(Comprehensive Sanction, size(3)) /// 
		   col(4) iscale(.8) xsize(12) ysize(7) imargin(-4 -4 -4 -4) 

gr export "Figures\Figure5_AllMargins_Q.eps", replace
gr export "Figures\Figure5_AllMargins_Q.pdf", replace 



